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We investigate the contribution of medium-modified jets and initial state fluctuations to the asym- 
metry in charged particle production with respect to the reaction plane. This asymmetry has been 
suggested as a compelling signature of the Chiral Magnetic Effect in QCD and make a study of 
"conventional" scenarios for the creation of such charged particle multiplicity fluctuations a timely 
endeavor. The different pathlength combinations of jets through the medium in non-central heavy 
ion collisions result in finite correlations of like and different charged particles emitted in the different 
hemispheres. Our calculation is based on combining jet events from YaJEM (Yet another Jet Energy 
Loss Model) and a bulk medium evolution. It is found that the jet production probabilities are too 
small to observe this effect. The influence of initial state fluctuations on this observable is explored 
using an event-by-event (3+1 )d hybrid approach that is based on UrQMD (Ultra-relativistic Quan- 
tum Molecular Dynamics) with an ideal hydrodynamic evolution. In this calculation momentum 
conservation and elliptic flow are explicitly taken into account. The asymmetries in the initial state 
are translated to a final state momentum asymmetry by the hydrodynamic flow profile. Dependent 
on the size of the initial state fluctuations the resulting charged particle asymmetries are in qual- 
itative agreement with the preliminary STAR results. The multi-particle correlation as proposed 
by the PHENIX collaboration can in principle be used to disentangle the different contributions, 
however is in practice substantially affected by the procedure to subtract trivial resonance decay 
contributions. 

PACS numbers: 25.75.-q,24.10.Lx,24.10.Nz,25.75.Gz 



I. INTRODUCTION 

Ultra-relativistic heavy-ion collisions at the Relativis- 
tic Heavy- Ion Collider (RHIC) are thought to have cre- 
ated a high temperature and pressure state of QCD mat- 
ter termed the Quark-Gluon-Plasma (QGP) 0^]: Two 
major discoveries stand out: (1) the emission of hadrons 
with a transverse momentum pt of several GeV/c or 
more is strongly suppressed (jet-quenching) , implying the 
presence of matter with a very large color opacity, and 
(2) the anisotropic ("elliptic") flow in non-central colli- 
sions is near the ideal hydrodynamic limit, requiring an 
early onset of the period during which the expansion is 
governed by fluid dynamics (earlier than 1 fm/c after the 
initial impact) as well as nearly ideal fluid properties with 
a viscosity-to-entropy density ratio rj/s <C 1. The mat- 
ter created at RHIC has been thus called the strongly 
interacting Quark-Gluon Plasma (sQGP) @ and quanti- 
fying its properties in terms of its transport coefficients 
is one of the major tasks of the current experimental and 
theoretical research program. 

One should note that QCD matter created in relativis- 
tic heavy-ion collisions may possess a very rich set of fea- 
tures, reflecting the fundamental symmetries (and viola- 
tions thereof) of QCD, well beyond the above noted char- 
acteristics of large opacity and near-ideal fluidity. Among 
the more intriguing features which are currently being 
searched for is the presence of the Chiral Magnetic Ef- 
fect (CME) [7HlQi: as was pointed out in Q , non-central 
heavy-ion collisions create a coherent magnetic field that 



may convert topological charge fluctuations in the QCD 
vacuum into global electric charge fluctuations with re- 
spect to the reaction plane. While the magnitude of the 
charge fluctuations which can be attributed to the CME 
is still a matter of ongoing debate HO]!, the STAR col- 
laboration has recently reported that charge fluctuations 
with respect to the reaction plane have indeed been seen 

;n]. 

In order to confirm the discovery of the CME in rel- 
ativistic heavy-ion collisions, one first needs to verify 
that experimental measurements cannot be understood 
in terms of " conventional" , well-established, phenomena 
observed at RHIC. It has recently been pointed out that 
correlations from charge and momentum conservation at 
freeze-out overlaid with elliptic flow may very well be 
able to describe the observed phenomena [IH, [l3| • In ad- 
dition, the experimental observable suggested to measure 
the CME has been argued to be predominantly sensitive 
to other effects [II 

In this paper we focus on estimating various conven- 
tional physics mechanisms that lead to charge parti- 
cle number asymmetry within a realistic model, namely 
medium-modified jets as well as initial state fluctuations, 
combined with well-known bulk physics such as momen- 
tum conservation, resonance decays and elliptic flow. In 
Section |H] the different observablcs that have been pro- 
posed are introduced. These observables are surely sen- 
sitive to fluctuations of multiplicity (and charge) if these 
fluctuations are correlated with the reaction plane. Both 
of the two szenarios that we are studying here, in-medium 
jets (in Section ITTT|) and initial state granularity (in Sec- 
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tion lIVp give a particular mechanism to generate multi- 
plicity fluctuations correlated in a peculiar way with the 
reaction plane. Only once all possible contributions to 
the observed charged particle fluctuations with respect 
to the reaction plane are accounted for, can we quantify 
the possible contribution of the CME to the measure- 
ments. 



II. DEFINITION OF OBSERVABLES 

The standard observable used to look for the CME is 
the angular correlation coefficient 
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where a and /3 represent positive and negative charges, 
(j) measures the angle with respect to the reaction plane 
orientation and M are the measured multiplicities. As ar- 
gued in [H, [TH, [TU , this observable is (contrary to initial 
expectations) dominated by a strong in-plane like-sign 
correlation which has no obvious relation with charge 
separation as proposed as signal of the CME. 

In the present study we consider two other newly devel- 
oped observables to quantify the correlations of charged 
particles with the reaction plane for which a computation 
using a detailed framework is feasible with a reasonable 
amount of CPU time. In addition, these observables seem 
to have a clearer relation to out-of-plane charge separa- 
tion. 
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and correspondingly for the LR-direction where N + 
and N~ are the numbers of positively and negatively 
charged particles in one event. 



The opposite charge correlation (A^ and ) are 
positive, if the particles tend to be aligned and negative 
for a charge separation scenario. Without any correlation 
this observable would be consistent with zero. The (A + ) 2 
and (A~) 2 are a measure of the overall charged particle 
fluctuation. The preliminary STAR results seem to be 
positive (on the order of 1 • 1CP 3 ) in the out-of-plane as 
well as the in-plane direction for mid-central Au+Au col- 
lisions at v /snn — 200A GeV. The main difference of this 
observable compared to the original angular correlation 
measurement [ll|, [l8[ is that each particle contributes 
with the same weight independent of its position in the 
respective hemisphere. 



B. Multi-Particle Correlation 

The PHENIX collaboration has recently proposed a 
different observable that is sensitive to a charge separa- 
tion with respect to the reaction plane as it is predicted 
from the chiral magnetic effect [l7| . This multi-particle 
correlation coefficient relies on an event- by- event mea- 
surement technique and is defined as follows 



A. Charge Asymmetry Observable 

The first observable we focus on for measuring the 
event-by-event charged particle asymmetry has been de- 
veloped by the Purdue group within the STAR collab- 
oration [16j. The measurement is based on counting 
charged particles that are emitted into different hemi- 
spheres with respect to the reaction plane. The differ- 
ent hemispheres are defined as follows: 'Up=U' means 
p y > 0, 'Down=D' refers to p y < 0, 'Left=L' stands for 
p x < and 'Right=R' for p x > 0. The UD-direction re- 
flects the out-of-plane direction whereas the LR-direction 
represents the in-plane direction. 

The charged particle asymmetry coefficients are de- 
fined in the following way 
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Here, the first sum runs over all positively charged par- 
ticles in one event whereas M is the number of negatively 
charged particles in one event and <fi is the azimuthal an- 
gle in the transverse plane in momentum space. A5 a ii is 
the analogous quantity where the ith particle is randomly 
chosen from the event. This observable is a multi-particle 
correlation that has to be measured event by event. The 
ratio of the charge sensitive distribution to the one for all 
charged particles is flat at 1, if there is no charge specific 
correlation (e.g. elliptic flow). A concave shape indi- 
cates charge separation out of plane (e.g. the presence 
of the CME) and a convex shape arises, if positively and 
negatively charged particles are preferably emitted in the 
same direction (e.g. resonance decays). 
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III. MEDIUM-MODIFIED JETS 
A. Overview 

The idea of jets as a mechanism to generate asymme- 
tries in the charged multiplicity distribution is based on 
the notion that in-medium energy loss from a leading par- 
ton is dominated by induced gluon radiation [TiJ |20| giv- 
ing rise to additional multiplicity production. Moreover, 
this mechanism is predominantly active out of plane: 
While radiative energy loss asymptotically has an L 2 de- 
pendence in a constant medium, the actual scaling of 
lost energy and produced multiplicity in a real situation 
is complicated and has no simple functional shape for 
a number of reasons, among them the spacetime evolu- 
tion of the bulk medium which reduces the density as a 
function of space and time and the re-interaction of radi- 
ated gluons with the medium, leading to a cascade of in- 
duced radiation till finite energy corrections become man- 
ifest However, the multiplicity of medium-induced 
hadron production is a monotonously rising function of 
both pathlength and medium density. 




FIG. 1: Sketch of the hadron multiplicity asymmetry induced 
by medium-modified showers: for a shower evolving in the 
medium, induced gluon radiation leads to larger multiplicity 
production. 

The basic idea how medium modified jets may gen- 
erate an asymmetry in the charged hadron multiplicity 
correlated out of the reaction plane is shown in Fig. [T] 
Let us start with a random initial hard vertex position in 
the transverse plane: if the vertex position is such that 
one jet is close to the surface and escapes unmodified 
whereas the other jet propagates a long distance through 
the medium, the average multiplicity of charged hadrons 
for both jets is expected to be very different. Further- 
more, this effect can be estimated to be stronger in the 
up-down direction than in the left-right direction, simply 
because the maximum variation in pathlength is larger 
for out of plane emission than for in-plane emission. 



Note that jets lead to an asymmetry in charged hadron 
production event by event even without a medium modi- 
fication, simply because the fragmentation of a hard par- 
ton to a hadron jet is a probabilistic process, and also be- 
cause jets in hadronic collisions may be due to a pQCD 
subprocess like qg — > qg in which the fragmentation of 
the gluon in the out state is much softer from the frag- 
mentation of the quark. However, in the presence of a 
medium, the medium-induced radiation is expected to 
dominate multiplicity production, and the vacuum case 
has no mechanism to account for a correlation with the 
reaction plane. In the following, we use the in-medium 
shower evolution code YaJEM (Yet another Jet Energy- 
loss Model) in the RAD (radiative energy loss) scenario to 
estimate the jet-induced event by event charged hadron 
multiplicity asymmetry. 

B. Model Description 

YaJEM is a Monte Carlo code to simulate the medium- 
modified evolution of a parton shower. It is based on 
the PYSHOW algorithm [22, 23] for shower evolution in 
vacuum, to which it reduces in the absence of a medium. 
A detailed description of the model is given in [3, HE| , 
here we just summarize the essential physics underlying 
the model. 

In YaJEM, the evolution from the initial parton to a 
final state parton shower is modelled as a series of branch- 
ing processes a — > b+c where a is called the parent parton 
and b and c are referred to as daughters. In QCD, the 
allowed branching processes are q — > qg, g — > gg and 
g — > qq. The kinematics of a branching is described in 
terms of the virtuality scale Q 2 and of the energy fraction 
z, where the energy of daughter b is given by E\, = zE a 
and of the daughter c by E c = (1 — z)E a . It is convenient 
to introduce t — InQ 2 / Aqcd where Aqcd is the scale 
parameter of QCD. t takes a role similar to a time in the 
evolution equations, as it describes the evolution from 
some high initial virtuality Qq (to) to a lower virtuality 
Qm (t m ) at which the next branching occurs. In terms 
of the two variables, the differential probability dP a for 
a parton a to branch is 

dP °- = Y\ ^Pa^bc(z)dtdz (7) 
b,c 

where the splitting kernels P a ^bc{z) read 

P^ 99 (z) = 4/3i±^ (8) 



P g ^ q - q {z) = N F /2(z 2 + (l-z) 2 ). (10) 

We do not consider electromagnetic branchings. Np 
counts the number of active quark flavours for given vir- 
tuality. The resulting system of equations describing the 
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branching processes in vacuum is solved numerically us- 
ing MC techniques utilizing the Pyshow code [H, [23[ . 

In order to make the link from momentum space where 
the shower evolution takes place to position space where 
the medium perturbations evolve, we assume that the 
average formation time of a shower parton with virtuality 
Q is developed on the timescale l/Q, i.e. the average 
lifetime of a virtual parton with virtuality Qb coming 
from a parent parton with virtuality Q a is in the rest 
frame of the original hard collision (the local rest frame 
of the medium may be different by a flow boost as the 
medium may not be static) given by 
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We assume that the actual formation time can then be 
obtained from a probability distribution 



P(r b ) = exp 
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which we sample to determine the actual formation time 
of the fluctuation in each branching. 

The medium modification in the RAD scenario appears 
as a change of parton virtuality AQ 2 dependent on a 
medium parameter q where 
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Here, the time r is given by Eq. (1121) 
is known in the simulation as the endpoint of the pre- 
vious branching process and the integration dC, is along 
the eikonal trajectory of the shower-initiating parton. If 
the parton is a gluon, the virtuality transfer from the 
medium is increased by the ratio of their Casimir color 
factors, 3/| = 2.25. q is a function of medium energy 
density e with the relation chosen such that the single 
pion suppression observed in 200 AGeV central Au-Au 
collisions at RHIC is described 
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As a result of the medium-induced change in virtuality, 
additional gluon radiation occurs which after hadroniza- 
tion leads to an increased multiplicity in the shower. 



C. Results 

In order to estimate the event by event charged hadron 
asymmetry, we embed back-to-back jet pairs of three dif- 
ferent energies (10, 20 and 40 GeV per jet) into a 3+1 
d hydrodynamical simulation (26j on a random position 
sampled from the binary collision profile. The jet pair 
then gets a random orientation with respect to the re- 
action plane. We then propagate both jets through the 
medium and evaluate Eq. (|13[) to determine the medium 
modification. After hadronizing the jets, we count the 
charged hadron content in the four different quadrants 
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FIG. 2: (Color online) Contribution of medium-modified jets 
to the charged particle asymmetry coefficients in the out-of- 
plane (left) and in-plane (right) direction for three different 
initial parton energies (10,20 and 40 GeV) for Au+Au colli- 
sions at -y/sNN = 200A GeV. Open symbols depict central col- 
lisions (b = 2.4 fm) while filled symbols represent mid-central 
collisions (b = 7.5 fm. 



and compute the asymmetry coefficients (A + ) 2 ,(A ) 2 
and . The results are shown in Fig. [2] 

The result indicates that jets are in principle capable 
of creating substantial charge correlations. In this sce- 
nario the charge correlation arises because of an asym- 
metry in the particle production. The correlation shown 
in Fig. [2] shows that in a developing shower the balancing 
charges are preferably emitted in the same hemisphere, so 
positively and negativley charged particles tend to align. 
However, somewhat surprisingly, no pronounced differ- 
ence between up-down and left-right is visible. This can 
be traced back to a number of reasons. First, Eq. (fT3|) 
evaluated in an expanding medium leads to small modi- 
fications at late times, i.e. for long paths, as q becomes 
very small when the energy density of the system is low. 
This means that in an expanding medium, partons do not 
actually probe the full length difference between in-plane 
and out-of-plane emission, but rather probe pathlength 
difference during the time period when modifications are 
significant only, i.e. the first ~ 3-5 fm/c. Given this 
'filter', the differences between in-plane and out-of-plane 
emission are considerably smaller than estimated naively 
based on geometry. Furthermore, especially for low px 
hadron production, the angular distribution of hadrons 
around the jet axis is wide [1^, [27j ■ Note that this can 
happen even if the jet shape itself is rather narrow, as 
the jet shape traces the flow of momentum and not mul- 
tiplicity. The fact that there can be significant multi- 
plicity created also in the down hemisphere even from 
a shower initiating parton travelling into the up hemi- 
sphere further weakens the difference between in-plane 
and out-of-plane emission. 
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So far, we have investigated the charge correlation co- 
efficients created by a single jet only. This is not a very 
realistic situation: First, high momentum jets of e.g. 40 
GeV are very rare events, so even if they create sizeable 
asymmetries in charged multiplicity, such events may not 
have much weight in an average over many events. Sec- 
ond, for low momentum jets, one cannot assume that 
only a single jet per event contributes. However, it is im- 
mediately clear that embedding multiple jets in an event 
weakens the asymmetry significantly. When we account 
for these effects by mixing jets with the probability of 
a jet with given energy being in the event determined 
from a LO pQCD calculation into a UrQMD background 
[H, [29| , the picture changes considerably: In Fig. [3] the 
contribution of jets appears as an insignificant correction 
to the background result. 
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FIG. 3: (Color online) Calculation of the charged parti- 
cle asymmetry coefficients in the out-of-plane (left) and in- 
plane (right) direction for two different centralities for Au+Au 
collisions at ^/ snn = 200A GeV. The open circles depict 
the results from the hadronic transport approach (UrQMD) 
that serves as a background for the full calculation including 
medium-modified jets that is represented by full squares. 

In hindsight, this could have been anticipated: While 
truly hard jets are too rare to influence the result, low 
Pt jets (below 5 GeV) tend to be created frequently, but 
the presence of multiple jets in an event tends to cancel 
their effect. Clearly, a bulk mechanism creating charge 
asymmetries is more efficient than a rare probe. 

Note also that it is apparent from Fig. [3] that UrQMD 
background results show that a pure hadronic transport 
approach is not capable of translating initial state fluctu- 
ations to the final state distributions effectively enough 
to lead to a significant signal as the A + /~ coefficients are 
on the order of 1 • 1CU 6 . 

In Fig. [4] the calculation of the multi-particle correla- 
tion coefficient is shown. Medium-modified jets, embed- 
ded in a UrQMD background lead to a convex shape of 
the ratio of the distributions in AS". The result for jets 
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FIG. 4: (Color online) Calculation of the multi-particle corre- 
lation coefficient in mid-central Au+Au collisions at ^/snn = 
2QQA GeV. The full squares depict the results from the full 
calculation including medium-modified jets to the hadronic 
transport approach (UrQMD) background event, whereas the 
full circles show the result only for particles with pr < 1 GeV. 

alone leads to the same shape as the combined events. 
But again, the probability for jet production is very low 
and therefore the results resemble very much the UrQMD 
results that are shown in Fig. |H1 Furthermore, the influ- 
ence of &pt < lGeV cut has been explored, since the pre- 
liminary measurements have been made in a limited px 
range and the chiral magnetic effect is mostly expected 
in the bulk low transverse momentum region. The low 
Pt result and the result for all charged particles are the 
same within error bars and especially the shape does not 
change. Since the distribution of the charge sensitive 
averaged sine differences has a higher peak around zero 
and lower flanks than the distribution that is not charge 
specific this result indicates that particles with opposite 
charges are preferably emitted in the same direction than 
in different directions. The underlying physics mecha- 
nism in both cases is correlated particle production from 
'clusters' - in one case jets, in the other decaying reso- 
nances. 



IV. INITIAL STATE FLUCTUATIONS 
A. Overview 

In the second part of this article the charged parti- 
cle asymmetry that arises from well-known bulk medium 
effects is investigated. The basic idea is that energy den- 
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sity clusters produced by the initial binary scatterings 
in heavy ion reactions might lead to asymmetry correla- 
tions in the particle production in different hemispheres. 
Due to the development of radial flow during the medium 
evolution the initial coordinate space granularity is trans- 
lated to a resulting momentum space anisotropy. These 
clusters are expected to enhance the multiplicity in a spe- 
cific hemisphere which results in alignment of the charged 
particle production instead of separation. In this picture, 
resonance decays might explicitly introduce a charge cor- 
relation of the outgoing particles, if e.g. a p meson decays 
into 7r + + 7r~ and there is enough collective flow, both 
of the pions will preferably be have a momentum in a 
similar direction. 

S. Pratt recently estimated the effects of elliptic flow, 
momentum and charge conservation on the observable 
charged particle correlation measurement [l3|. Our ap- 
proach follows a similar spirit but the effect of initial 
state fluctuations is emphasized in addition. This kind 
of event-by-event fluctuations has been found to be cru- 
cial also for other two-particle correlation measurements 
(the 'ridge') at the highest RHIC energies [30j . 



where N = (^r) 2 -pr^ — E c [ provides the proper nor- 
malisation, e c f and E c { are the energy density and to- 
tal energy of the particle in the computational frame, 
while (xp,y p ,Zp) is the position vector of the particle. 
Summing over all single particle distribution functions 
leads to distributions of energy-, momentum- and baryon 
number-densities in each cell. 

The transition from the hydrodynamic evolution to the 
transport approach when the matter is diluted in the 
late stage is treated as a gradual transition on an ap- 
proximated iso-eigentime hyper-surface (see [37], [38| for 
details). The final rescatterings and resonance decays are 
taken into account in the hadronic cascade. 

To investigate the effect of different initial state gran- 
ularities the width of the Gaussian distribution (fT4)) is 
varied from CQauss = 0.8 — 2.0 fm. Please note that this 
is still an event- by- event setup and in this respect differ- 
ent from averging over many fluctuating initial conditions 
as has been done in [39$ . In a hydrodynamic calculation 
starting from smooth initial profiles the correlation ob- 
servable would be zero, since there are no event- by- event 
fluctuations included. 



B. Model Description 



The charged particle asymmetry that is generated 
by clusters in the initial energy density distribution is 
quantified by using a hybrid approach based on the 
UrQMD hadronic transport approach including a (3+1)- 
dimensional one fluid ideal hydrodynamic evolution [3lL 
[32l | for the hot and dense stage of the reaction while the 
early non-equilibrium stage and the final decoupling is 
treated in the hadronic cascade [HI, [34]] [I(| • 

The event-by-event-setup that is naturally employed 
in this approach is crucial for the observation of charged 
particle asymmetries. For Au+Au collisions at the high- 
est RHIC energy the starting time for the hydrodynamic 
evolution has been chosen to be t s tart = 0.3 — 1.8 fm 
depending on the amount of fluctuations to fit the final 
pion multiplicity at midrapidity, since this is crucial for 
the present study and allows us to investigate the effect 
of the initial state granularity on the charge correlations. 
How the change in the starting time affects other bulk ob- 
servables like spectra and elliptic flow will be discussed 
in a further publication |35|. Only the matter around 
midrapidity (\y\ < 2) is considered to be locally thermal- 
ized and takes part in the ideal hydrodynamic evolution. 
The more dilute spectator/corona regions are treated in 
the hadronic cascade approach. To map the point parti- 
cles from the UrQMD initial state to energy, momentum 
and net baryon density distributions three-dimensional 
Gaussian distributions are used that represent one parti- 
cle each 13611. 




e cl (x,y,z) = TVexp- 



(x - x p ) 2 + (y- y v ) 2 + (7z(z - z v )) 2 
2a 2 
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FIG. 5: (Color online) Energy density distribution in the 
transverse plane for a central (b=0 fm) Au+Au collision at 
^/snn = 200^4 GeV with <TGauss = 0.8 fm and tstart = 0.3 fm. 

Fig. [5][7] show the initial energy density distribution 
as it is obtained from UrQMD for three different val- 
ues of the Gaussian width ccauss = 0.8, 1.0 and 2.0 fm. 
A smaller width leads to entropy production during the 
ideal hydrodynamic evolution indicating that the solu- 
tion of the differential equation is not stable anymore. 
Therefore, the highest fluctuations are generated with the 
smallest width of 0.8 fm, while for larger widths the ini- 
tial fluctuations are more and more smeared out. Smaller 
characteristic sizes could be physically motivated by a 
color flux tube picture where structures in the transverse 
plane on the order of 0.3 fm would be expected. Our 
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parameter choices stay in a regime where the implemen- 
tation is technically stable and the characteristic sizes 
correspond to typical sizes of a proton. 



b=0 fm; |y|<2.0; t start =0.5 fm; G=1.0 fm 
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on event-by-event particle production in different hemi- 
spheres. Initial state clusters and a coordinate space 
-momentum space correlation due to radial flow leads 
to observable differences in the charged particle pro- 
duction in different parts of the system (see Fig. [8]). 
The UD and LR charged particle asymmetry coefficients 
have been calculated for four different initial granulari- 
ties in the hybrid approach. The open symbols depict 
the results for central (b < 3.4 fm) Au+Au collisions 
at ^/snn = 200A GeV and the filled symbols represent 
mid-central (b = 5 — 9 fm) Au+Au collisions. The overall 
fluctuations of charged particles are larger in non-central 
collisions than in central collisions because the overall 
multiplicity is smaller. The opposite charge asymmetry 
is positive and similar in the out-of-plane and the in- 
plane direction. This confirms the expectations that in 
this scenario the particles tend to align independent of 
their charge. Within our limited statistics, the absolute 
value of the difference between UD and LR can be con- 
strained to be smaller than f • 10 -3 . 



FIG. 6: (Color online) Energy density distribution in the 
transverse plane for a central (b=0 fm) Au+Au collision at 
a/snn = 200A GeV with crGauss = 1.0 fm and tstart = 0.5 fm. 
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FIG. 7: (Color online) Energy density distribution in the 
transverse plane for a central (b=0 fm) Au+Au collision at 
a/snn = 200A GeV with aGauss = 2.0 fm and tstart = 1.8 fm. 



C. Results 

The hybrid approach provides the opportunity to 
study the influence of initial state fluctuations that can 
be varied in size without adjusting any other parameter 
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FIG. 8: (Color online) Calculation of the charged particle 
asymmetry coefficients in the out-of-plane (left) and in-plane 
(right) direction for four different values of crGauss. The open 
symbols depict the results for central (6 < 3.4 fm) Au+Au col- 
lisions at y / snn = 200A GeV and the filled symbols represent 
mid-central (6 = 5 — 9 fm) Au+Au collisions. 

Higher granularity in the initial state leads to higher 
opposite charge asymmetry coefficients that are on the 
order of the preliminary STAR data 16]. This physics 
mechanism will contribute to the background of a po- 
tential measurement of the chiral magnetic effect and 
needs to be taken into account, before trying to look for 
a potentially much smaller signal. On the other hand, 
the comparison of hybrid calculation with different ini- 
tial state granularities to event- by- event correlation and 
fluctuation measurements might be helpful to constrain 
the structure of the initial state of relativistic heavy ion 
reactions. 

Fig. [S] shows the result for the multi-particle correla- 
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FIG. 9: (Color online) Calculation of the multi-particle corre- 
lation coefficient in mid-central Au+Au collisions at ^/snn = 
200A GeV. The open symbols depict the results from the 
hadronic transport approach (UrQMD) (circles indicate an 
additional pr < 1 GeV cut) and the filled symbols represent 
calculations within the hybrid approach. 



tion coefficient that has been defined in Section lHBI from 
the UrQMD hadronic transport approach (open symbols) 
and the hybrid approach (filled symbols). It has been 
checked that the results for this observable do not depend 
on the choice of cxGauss- The low transverse momentum 
region shows the same behaviour as without px cut. If 
one reduces the number of particles by acceptance (or 
other) cuts, the distribution broadens, however does not 
change qualitatively in shape. Again the correlations in- 
duced by these approaches lead to a convex shape and 
can be attributed to resonance decays whereas a charge 
separation scenario as expected from the chiral magnetic 
effect leads to a concave shape in this observable. The 
initial state fluctuations have to be considered as a non- 
trivial source of correlations of the particles emitted in a 
heavy ion reaction. 

It is likely that any measurement of the observable is 
dominated by resonance decay effects which are known 
to be present in nature. In order to observe a different 
mechanism (like the CME) at work, this distribution has 
to be subtracted. In order to get an idea of the sen- 
sitivity of the results to the subtraction, we have tried 



to isolate the effect of initial state fluctuations by sub- 
tracting the UrQMD distributions in AS from the hybrid 
model distributions. It turns out that the procedure is 
very sensitive to the subtraction - we were able to ob- 
serve even qualitative differences in shape which could 
potentially mock up the CME signal. 

V. SUMMARY AND CONCLUSIONS 

We have investigated charge asymmetry and multi- 
particle correlations with the reaction plane generated in 
rclativistic heavy-ion collisions at RHIC within a realistic 
transport approach containing a full medium evolution 
and medium-modified jets. While jets are in principle 
capable of generating the aforementioned multi-particle 
correlations, we were able to rule them out as a signifi- 
cant contributor to such correlations due to the low jet 
production rate. However, we find that initial state fluc- 
tuations can generate values for the charge asymmetry 
coefficients on the same order of magnitude as seen by 
experiment. Our results were obtained within a hybrid 
transport approach, based on an event-by-event hydrody- 
namic evolution, with fluctuating initial conditions and 
microscopic hadronic afterburner for the realistic treat- 
ment of hadronic final state effects, including resonance 
decays. The precision of our calculation was not suf- 
ficient to quantify the difference between in-plane and 
out-of-plane, due to the very large computational effort 
involved. However, the coefficients themselves show no 
sign of charge separation. The mult i- particle correlation 
measurement proposed by the PHENIX collaboration ap- 
pears to be potentially the most suitable of the studied 
observables for differentiating a CME from other mech- 
anisms. However, due to the need of subtracting known 
background effects like resonance decays, it is currently 
unclear whether the predicted concave shape as a signal 
for the CME can be robustly established. 



VI. ACKNOWLEDGEMENTS 

We are grateful to the Open Science Grid and the Cen- 
ter for the Scientific Computing (CSC) at Frankfurt for 
the computing resources. The authors thank Fuqiang 
Wang, Roy Lacey and Berndt Miiiler for helpful discus- 
sions and Dirk Rischke for providing the f fluid hydro- 
dynamics code. H.P. acknowledges a Feodor Lynen fel- 
lowship of the Alexander von Humboldt foundation. This 
work was supported in part by U.S. department of Energy 
grant DE-FG02-05ER41367 and by project 133005 of the 
Academy of Finland. T.R. is supported by an Academy 
Research Fellowship of the Academy of Finland. 



[1] Arsene I et al. (BRAHMS) 2005 Nucl. Phys. A757 1-27 283 

[2] Adcox K et al. (PHENIX) 2005 Nucl. Phys. A757 184- [3] Back B B et al. 2005 Nucl. Phys. A757 28-101 



9 



Adams J et al (STAR) 2005 Nucl. Phys. A757 102-183 
Muller B and Nagle J L 2006 Ann. Rev. Nucl. Part. Set. [24 
56 93-135 [25 
Gyulassy M and McLerran L 2005 Nucl. Phys. A750 30- [26 
63 

D. E. Kharzeev, L. D. McLerran and H. J. Warringa, [27 
Nucl. Phys. A 803, 227 (2008) [28 
K. Fukushima, D. E. Kharzeev and H. J. Warringa, Phys. 
Rev. D 78, 074033 (2008) [29 
K. Fukushima, D. E. Kharzeev and H. J. Warringa, Phys. [30 
Rev. Lett. 104, 212001 (2010) 

M. Asakawa, A. Majumder and B. Muller, Phys. Rev. C 

81, 064912 (2010) ' [31 

B. I. Abelev et al. [STAR Collaboration], Phys. Rev. Lett. 

103, 251601 (2009) [32 

S. Schlic hting and S. Prat t, larXiv:1005. 53411 [nucl-th], 
S. Pratt. laFXiv:1002.1758l [nucl-th]. [33 
F. Wang, Phys. Rev. C 81, 064902 (2010) 

A. Bzdak, V. Koch and J. Liao, Phys. Rev. C 81, 031901 [34 
(2010) 

F. Wang, talk given at 'P- and CP- odd Effects in Hot [35 
and Dense Matter' on behalf of the STAR collaboration, 
April 2010. [36 
N. Ajitanand, talk given at the annual AGS-RHIC user 
meeting 2010 on behalf of the PHENIX collaboration. 

B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C [37 
81, 054908 (2010) 

R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and [38 
D. Schiff, Nucl. Phys. B 484, 265 (1997) 
B. G. Zakharov, JETP Lett. 65, 615 (1997) 
R. B. Neufeld and T. Renk, Phys. Rev. C 82, 044903 [39 
(2010) 

M. Bengtsson and T. Sjostrand, Phys. Lett. B 185, 435 [40 
(1987). 

M. Bengtsson and T. Sjostrand, Nucl. Phys. B 289, 810 



(1987). 

T. Renk, Phys. Rev. C 78, 034908 (2008) 
T. Renk, Phys. Rev. C 79, 054906 (2009) 

C. Nonaka and S. A. Bass, Phys. Rev. C 75, 014902 
(2007) 

T. Renk, Phys. Rev. C 80, 044904 (2009) 

S. A. Bass et al, Prog. Part. Nucl. Phys. 41, 255 (1998) 

[Prog. Part. Nucl. Phys. 41, 225 (1998)] 

M. Bleicher et al, J. Phys. G 25, 1859 (1999) 

J. Takahashi, B. M. Tavares, W. L. Qian, F. Grassi, 

Y. Hama, T. Kodama and N. Xu, Phys. Rev. Lett. 103, 

242301 (2009) 

D. H. Rischke, S. Bernard and J. A. Maruhn, Nucl. Phys. 
A 595, 346 (1995) 

D. H. Rischke, Y. Pursun and J. A. Maruhn, Nucl. Phys. 

A 595, 383 (1995) [Erratum-ibid. A 596, 717 (1996)] 

H. Petersen, J. Steinheimer, G. Burau, M. Bleicher and 

H. Stocker, Phys. Rev. C 78, 044901 (2008) 

H. Petersen and M. Bleicher, Phys. Rev. C 79, 054904 

(2009) 

H. Petersen, C. Coleman-Smith, S. A. Bass and 

R. Wolpert, larXiv: 1012.46291 [nucl-th] . 

J. Steinheimer, M. Bleicher, H. Petersen, S. Schramm, 

H. Stocker and D. Zschiesche, Phys. Rev. C 77, 034901 

(2008) 

Q. f. Li, J. Steinheimer, H. Petersen, M. Bleicher and 

H. Stocker, Phys. Lett. B 674, 111 (2009) 

J. Steinheimer, V. Dexheimer, H. Petersen, M. Bleicher, 

S. Schramm and H. Stoecker, Phys. Rev. C 81, 044913 

(2010) 

H. Petersen and M. Bleicher, Phys. Rev. C 81, 044906 
(2010) 

The code is available as UrQMD-3.3pl at 
http://urqmd. org] 



